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FINAL  REPORT 

PLASMA  WAVE  ELECTRONICS 
MICHAEL  S.  SHIJR 

Electrical,  Computer  and  Systems  Engineering  Department 
Rensselaer  Polytechnic  Institute,  Troy,  NY  12180-3590,  USA 


Plasma  waves  are  oscillations  of  electron  density  in  time  and  space.  In  deep  submicron  field  effect 
transistors  plasma  wave  frequencies  lie  in  the  terahertz  range  and  can  be  tuned  by  applied  gate  bias.  Since 
the  plasma  wave  frequency  is  much  larger  that  the  inverse  electron  transit  time  in  the  device,  it  is  easier  to 
reach  “ballistic”  regimes  for  plasma  waves  than  for  electrons  moving  with  drift  velocities.  In  the  ballistic 
regime,  no  collisions  of  electrons  with  impurities  or  lattice  vibrations  occur  on  a  time  scale  on  the  order  of 
the  plasma  oscillation  period,  and  the  device  channel  acts  as  a  resonant  cavity  for  the  plasma  waves, 
making  possible  tunable  resonant  detection  or  even  emission  of  the  electromagnetic  radiation  in  the 
terahertz  range.  In  this  report,  we  present  the  theory  of  plasma  waves  in  field  effect  transistors;  discuss 
instabilities  of  these  waves  in  different  device  structures  and  their  applications  for  detection  and 
generation  of  the  terahertz  radiation.  The  summary  of  this  report  will  be  published  in  Reference 


1. Introduction 

The  terahertz  range  of  frequencies  is  often  referred  to  as  the  “terahertz  gap”,  since  this  frequency  range  lies  in 
between  the  frequency  ranges  of  electronic  and  photonic  devices.  The  upper  frequency  limit  of  transistors 
operating  in  conventional  regimes  is  limited  by  the  transit  time  of  carriers  under  the  gate  (for  a  field  effect 
transistor)  or  across  the  base  and  collector  depletion  region  (for  a  bipolar  junction  transistor).  Scaling  of  feature 
sizes  of  silicon  Metal  Oxide  Semiconductor  Field  Effect  Transistors  (MOSFETs)  and  compound  semiconductor 
Heterostructure  Field  Effect  Transistors  (HFETs)  and  Heterojunction  Bipolar  Transistors  has  pushed  the  device 
parameters  into  the  region,  where  the  transistor  operation  at  a  few  hundred  gigahertz  became  feasible.  1 
However,  device  feature  sizes  approach  the  values,  such  that  fundamental  physics  limitations  lead  to 
diminishing  returns  on  investment  in  further  scaling,  and  the  transit  time  limited  regimes  face  these  limitations 
in  trying  to  approach  the  terahertz  range  of  frequencies.  On  the  other  hand,  since  the  quanta  of  terahertz 
radiation  are  much  smaller  than  the  thermal  energy  at  room  or  even  at  liquid  nitrogen  temperatures,  photonic 
devices  using  interband  or  intersubband  transitions  have  to  operate  at  cryogenic  temperatures  (see  Fig.  1). 

Plasma  waves  are  oscillations  of  electron  density  in  time  and  space,  and  in  deep  submicron  field  effect 
transistors,  typical  plasma  frequencies,  o)P,  lie  in  the  terahertz  range  and  do  not  involve  any  quantum  transitions. 
Hence,  using  plasma  wave  excitation  for  detection  and/or  generation  of  terahertz  oscillations  is  a  very 
promising  approach,  and,  as  illustrated  by  Fig.  1,  the  terahertz  gap  in  the  electromagnetic  spectrum  between 
electronic  and  photonic  devices  can  be  closed  by  plasma  wave  electronics  devices.  We  call  this  approach 
plasma  wave  electronics, 2A 


1 M.  S.  Shur  and  V.  Ryzhii,  Plasma  wave  electronics,  in  Terahertz  Sensing  Technology,  Volume  I.  Electronic  Devices  &  Advanced 
Technology,  Dwight  Woolard,  William  Loerop,  mid  Michael  S.  Shur,  Editors,  World  Scientific  (2003),  in  press 


The  properties  of  plasma  waves  depend  on  the  electron  density  and  on  the  dimensions  and  geometry  of  the 
electronic  system.  In  a  three  dimensional  case,  the  plasma  oscillation  frequency  is  nearly  independent  of  the 
wave  length,  in  a  gated  two-dimensional  electron  gas  (2DEG),  the  plasma  wave  have  a  linear  dispersion  law 
similar  to  that  of  sound  waves  or  light  in  vacuum.  In  this  case,  the  plasma  wave  velocity,  s,  is  proportional  to 
the  square  root  of  the  electron  sheet  density  and  can  be  easily  tuned  by  the  gate  bias  that  controls  the  2DEG 
density. 
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Fig.  1 .  Terahertz  range 


This  is  the  physical^  basis  for  plasma  wave  electronics  devices  that  involve  two  dimensional  electrons  in  a 
field  effect  transistor  or  in  quantum  well  structures. 5,6  In  deep  submicron  devices,  the  electron  transit  time 
might  become  comparable  to  or  even  smaller  than  the  electron  momentum  relaxation  time,  xm.  Under  such 
conditions,  the  electron  transport  approaches  the  ballistic  regime.7  (Recently,  the  Lucent  group  reported  on  a 
ballistic  transistor ,  and  Kawaura  et  al  demonstrated  a  14  nm  transistor.^)  For  electrons  moving  in  a  device 
channel  with  a  certain  drift  velocity,  there  are  two  ramifications  of  the  ballistic  mode  of  transport. 

In  low  electric  fields,  the  effective  electron  mobility  in  short  channel  (submicron)  devices  might  be  much 
icm*m2  r  ^an  ^*e  e'ec^ron  m°bil ity  in  long  channel  devices.  This  reduction  was  predicted  for  ballistic  devices  in 
,  ,  ,  and,  as  was  shown  recently13,  is  consistent  with  experimental  data  for  submicron  AlGaAs/GaAs  High 
Electron  Mobility  Transistors  (HEMTs).  This  ballistic  effect  should  be  also  important  in  short  channel  devices 
implemented  in  silicon,  silicon-germanium,  AlN/GaN/InN  or  any  other  semiconductors.  The  physical  reasons 
for  a  drastic  mobility  reduction  are  related  to  the  ballistic  motion  first  predicted  in  1979. 6  In  ballistic  field  effect 
transistors,  electrons  travel  from  the  source  to  the  drain  without  any  collisions  with  impurities  or  phonons,  and 
electrons  propagate  in  the  device  channel  with  a  randomly  oriented  thermal  velocity. 


v 


th 


rUBT' 

mt  j 


(1) 


(where  m  is  the  effective  mass,  kg  is  the  Boltzmann  constant,  and  T  is  temperature)  or  with  a  Fermi  velocity,  vF, 
for  a  degenerate  electron  gas.  Hence,  electrons  have  only  a  limited  time  to  accelerate  in  the  electric  field  and 
acquire  a  drift  velocity.  This  acceleration  time  is  determined  by  L/vth,  where  L  is  the  device  length,  (or  by  L/vF 
in  a  degenerate  case).  As  a  result,  in  low  electric  fields,  the  current  is  proportional  to  the  electric  field  and  to  the 
electron  concentration,  just  like  in  the  collision-dominated  case,  but  the  electron  mobility  has  to  be  substituted 
by  an  effective  "ballistic’  mobility,  which  (for  a  non-degenerate  case)  is  given  by  1 1 

2  qL 

M ballistic  —  (2) 
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In  high  electric  fields,  electron  velocity  in  a  ballistic  device  is  expected  to  be  higher  or  even  much  higher 
that  in  a  collision  dominated  device  7,  and  the  bulk  plasma  oscillations  should  lead  to  space  oscillations  of  the 
electron  density  and  resulting  S-type  current-voltage  characteristic. 14 


Plasma  oscillations  in  a  field  effect  transistor  are  affected  by  streamlines  of  the  electric  field  directed  from 
the  channel  toward  the  gate,  and  the  velocity  of  plasma  waves  in  a  field  effect  transistor  is  typically  much 
higher  (by  an  order  of  magnitude  or  more)  than  the  electron  drift  velocity.  Hence,  the  characteristic  transit  time 
is  much  shorter,  and  the  condition  0)pXm  >  1  is  much  easier  to  meet  that  the  condition  2nxm/tIr>  1,  where  t,r  is  the 
electron  transit  time.  When  the  condition  (OpXm  >  1  is  met,  a  channel  of  a  field  effect  transistor  can  serve  a 
resonant  cavity  for  the  plasma  waves.  As  discussed  below,  the  fundamental  frequency  of  this  cavity  can  be 
tuned  by  changing  the  gate  bias,  and  a  high  mobility  deep  submicron  FET  can  be  used  for  resonant  detection, 
mixing,  multiplication,  5  and  even  generation  16  of  terahertz  radiation. 

We  will  start  from  a  general  discussion  of  plasma  waves  for  different  geometries  with  a  special  emphasis  on 
the  plasma  waves  in  a  FET  channel  used  for  plasma  wave  electronics.  We  will  then  consider  the  analogy 
between  plasma  waves  in  gated  two-dimensional  electron  gas  (2DEG)  and  shallow  water  waves  and  the  analogy 
between  plasma  waves  in  ungated  2DEG  and  deep  water  waves.  This  will  be  followed  by  the  analysis  of  the 
instability  of  plasma  waves  in  gated  2DEG  channels.  We  will  then  review  experimental  results  on  resonant  and 
non-resonant  detection  and  new  ideas  on  using  resonant  tunneling  structures  for  sharpening  surface  plasmon 
resonances. 


2.Plasma  Wave  Oscillations. 


Following  17,  let  us  consider  plasma  wave  oscillations  in  the  systems  of  different  dimensions  by  neglecting 
collisions  and  considering  only  the  average  drift  velocity,  v.  The  dispersion  relations  for  plasma  waves  can  be 
derived  from  the  small  signal  equation  of  motion  and  the  continuity  equation,  which,  under  these  assumptions, 
are  given  by: 

0j  e2n 

—  =  E — .  m 


—  +  div\  =  0 .  (4) 

where  j  =  qn\  is  the  current  density,  e  is  the  electronic  charge,  n  is  the  electron  density,  m  is  the  electronic 
effective  mass,  and  p  is  a  small-signal  charge  density  related  to  a  deviation  of  n  from  its  equilibrium  value,  E  is 
the  small  signal  electric  field.  Eq.  (3)  follows  from  the  Newton  Second  Law  of  Motion,  where  electron 
scattering  is  neglected.  For  a  3D  geometry,  the  j,  w,  and  p  are  current  per  unit  area,  electron  concentration,  and 
electric  charge  per  unit  volume,  respectively,  whereas  in  the  2D  case,  the  j,  n,  and  p  are  ,  whereas  in  the  2D 

case,  the  j,  n,  and  p  are  current  per  unit  length,  electron  concentration  per  unit  area,  and  the  electric  charge  per 
unit  area,  respectively. 

Differentiating  Eq.  (4)  with  respect  to  time  and  using  Eq.  (3),  we  obtain: 

d2p  q2n  . 

_f  +  i_rf,vE  =  0.  (5) 

For  the  three  dimensional  case,  the  relation  between  E  and  p  is  given  by: 

divE  =  — .  (6) 

£ 

where  e  is  the  dielectric  permittivity.  Substituting  Eq.  (6)  into  Eq.  (5),  we  obtain  the  well-known  expression  for 
the  3D  plasma  frequency 


where  s  is  the  dielectric  permittivity  of  a  semiconductor. 

For  electrons  under  the  gate  of  an  FET,  the  electron  sheet  concentration  in  the  channel  is  proportional  to  the 
voltage  difference  between  the  gate  and  channel  potentials: 

qn  =  CU  .  (8) 

Here  U  =  Ug  -  Uc  -  Ut,  Ut  is  the  threshold  voltage,  Ug  -  Uc  is  the  potential  difference  between  the  gate  and  the 
channel,  C  =  e/d  is  the  gate-to-channel  capacitance  per  unit  area,  and  d  is  the  gate-to-channel  separation.  This 


equation  is  valid  when  U  changes  along  the  channel  on  the  scale  large  compared  to  d  (so-called  gradual  channel 
approximation).  From  Eq.  (8),  we  find  the  relationship  between  the  electric  field  E  =  -V  U  and  two- 
dimensional  charge  density  p  =  q(n  -  n0) 


E  =  --Vp.  (9) 

Substituting  Eq.  (7)  into  Eq.  (3),  we  obtain  the  two-dimensional  wave  equation  for  the  surface  charge  p 

d2  p  2 

tfsAP  =  °-  00) 

where  A  is  the  two-dimensional  Laplace  operator,  and 

s=M=M-  <u) 

is  the  velocity  of  the  surface  plasma  waves  that  have  a  linear  dispersion  law: 

np=sk.  (12) 

Here  k  is  the  wave  vector  of  the  plasma  waves.  Using  Eq.  (6)  and  (7),  we  can  express  the  plasma  wave 
velocity  in  terms  of  the  gate  voltage  swing: 

In  a  similar  way,  using  Eq.  (6)  and  the  relation  between  the  in-plane  electric  field  and  surface  charge 
density,  one  can  obtain  the  dispersion  law  for  the  plasma  waves  for  ungated  2D  electron  gas 


and  for  the  plasma  waves  propagating  along  a  one-dimensional  wire: 


con  =  s.k  In — 

P  1 


Here  r  is  the  radius  of  the  wire,  s\  is  the  velocity  of  the  one-dimensional  plasma  waves  given  by 


nq 

4mt£ 


The  case  of  a  gated  wire  is  very  similar  to  that  of  a  FET  channel: 


<°P  ~  su ^ 


where 


q2n_  fqU 


C\d  =  .  (18) 

The  equations  for  the  ungated  2D  gas  are  valid  when  kd  «  1  (in  the  opposite  limiting  case  kd»  1  the  existence 
of  the  metallic  gate  is  irrelevant,  and  one  obtains  the  dispersion  relation  given  by  Eq.  (13)).  For  the  ID  case  Eq. 
(16)  holds  when  kr  «  1 . 

Fig.  2  gives  the  summary  of  the  plasma  wave  frequencies  and  dispersion  relations  for  systems  of  different 
dimensions  and  geometry. 


3D 


Fig.  2.  Plasma  wave  frequencies  for  different  sample  geometries. 3 
Fig.  3  shows  the  dispersion  relations  for  different  geometries.  (The  dispersion  curve  calculations  for  Fig.  3  used 
the  following  parameters:  dielectric  constant,  sr  =  12.9,  effective  mass,  mr  =  0.067  m0;  3D  electron 
concentration,  —  5x10  m  ,  2D  electron  concentration,  «2d  =  1.5x1 012  m'2,  ID  electron  concentration,  = 
1 .41x10  m  ,  gate-to-channel  separation  for  the  1 D  gated  structures,  d0~  3  nm,  quantum  wire  radius,  r  =  3  nm.) 


Fig.  3.  Plasma  wave  dispersion  relations  for  different  geometries. 

The  similarity  between  the  dispersion  relations  for  the  gated  2DEG  case  and  gated  and  ungated  1DEG  case 

implies  that  the  results  discussed  below  for  the  gated  2DEG  should  equally  apply  to  gated  and  ungated  quantum 
wires. 

In  ultra  quantum  limit  in  quantizing  magnetic  fields,  it  is  possible  to  achieve  a  situation  when  electrons 
propagate  as  ID  particles  but  electrostatic  interaction  is  three  dimensional.  This  happens  if  the  electron  energy 
is  less  than  the  separation  between  the  Landau  subbands  and  optical  phonon  energy.  In  such  a  system,  quasi¬ 
elastic  collisions  of  electrons  (with  impurities  and  acoustic  phonons)  can  only  invert  the  direction  of  the  electron 
propagation.  This  can  minimize  the  effect  of  these  collisions  on  plasma  instabilities.  As  an  example,  Ivanov  and 
Ryzhii  considered  the  instability  of  the  electron  plasma  optically  generated  in  the  lowest  Landau  subband.18  The 
dispersion  equation  for  this  system  coincides  with  that  for  two  colliding  electron  beams  both  in  a  classical  case 
(for  small  wave  numbers)  and  in  a  quantum  case  (when  Wigner  equation  should  be  used).  It  is  instructive  that 
elastic  scattering  does  not  influence  the  instability  threshold. 


The  same  authors  (jointly  with  Yu.  Sigov  and  V.  Kustov)  modeled  the  dynamics  of  this  instability  (which 
they  called  the  collective  relaxation)  using  a  ensemble  particle  method  and  observed  fast  transformation  of  the 
distribution  function  due  to  interaction  with  the  self-consistent  electric  field  (unpublished).  This  instability  led 
to  a  strong  turbulence  since  its  increment  was  close  to  the  instability  frequency. 

3.  Plasma  Waves  in  FETs 

Chaplik  .  was  the  first  to  consider  plasma  waves  in  a  gated  2D  electron  gas  for  electrons  on  the  surface  of 
liquid  helium.  Nakayama.  analyzed  plasma  waves  for  a  FET.  Allen  et  al.  1  observed  infrared  absorption  and 
Tsui  et  al  observed  weak  infrared  emission  related  to  such  waves  in  silicon  inversion  layers.  More  recent 
studies  by  Burke  et  al  J  of  high  mobility  AlGaAs/GaAs  gated  heterostructures  revealed  the  resonance 
impedance  peaks  related  to  the  plasma  waves.  Dyakonov  and  Shur  used  hydrodynamic  equations  to  analyze 
plasma  waves  in  2DEG  and  predicted  the  instability  of  plasma  waves  in  a  high  mobility  field  effect  transistor. 16 
They  also  developed  the  theory  of  plasma  wave  electronics  devices,  such  as  detectors  and  mixers.  15  Lu  et  al 
reported  the  evidence  of  the  resonant  plasma  wave  detection  by  a  FET  at  the  third  harmonic. 24  Recently,  Knap 
et  al  reported  on  the  resonant  plasma  wave  detection  by  a  FET  at  the  fundamental  harmonic,  and  Peralta  et  al 
reported  on  a  similar  detection  in  multi  gated  periodic  structures  with  2DEG,26 

The  gate  electrode  in  a  FET  (see  Fig.  4)  is  separated  from  the  channel  by  the  gate  insulator  (which  is  a 
doped  or  undoped  wide  band  gap  semiconductor,  such  as  AlGaAs  in  a  typical  High  Electron  Mobility 
Transistors  and  Si02  in  silicon  MOSFETs). 


Fig.  4.  Field  Effect  Transistor 


A  FET  operates  in  two  different  regimes.  At  gate  biases,  VG,  smaller  than  the  FET  threshold  voltage,  VT,  the 
source  and  drain  contacts  are  separated  by  a  potential  barrier  and  the  current  between  the  source  and  the  drain  is 
very  low.  An  increase  in  the  gate  voltage  decreases  the  height  of  this  barrier  leading  to  an  exponential  rise  in 
current.  When  VG  >VT,  electrons  attracted  by  a  more  positive  gate  charge  form  a  narrow  sheet  of  mobile  charge 
at  the  interface  between  the  semiconductor  and  the  gate  insulator.  This  electron  sheet  and  the  gafo  contact  form 
a  capacitor,  which  capacitance  per  unit  area  is  where  C  =  s/d,  where  d  is  the  gate-to  channel  separation,  Above 
the  threshold,  the  surface  concentration,  ns  in  the  FET  channel  is  given  by  (see  Eq.  (E) 

ns  =  CU  /  q  (19) 

As  mentioned  above,  Eq.  (19)  represents  the  usual  gradual  channel  approximation,  which  is  valid  when  the 
characteristic  scale  of  the  potential  variation  in  the  channel  is  much  greater  than  the  gate-to-channel  separation. 


The  equation  of  motion  (the  Euler  equation)  is 


dv  dv  q  dU  v 

— +v— +— — — + —  =  0, 
dt  dx  m  8x  r 


(20) 


where  dU/dx  is  the  longitudinal  electric  field  in  the  channel,  v(x,t)  is  the  local  electron  velocity,  and  xm  is  the 
momentum  relaxation  time.  The  usual  continuity  equation  can  be  written  as: 

t  akv)_0 

dt  dx 


(21) 


and  re-written  as 


mj  rtuv) 

T+  *  0  <22> 

taking  into  account  Eq.  (19). 

In  the  limiting  case  of  a  ballistic  FET,  xm  tends  to  infinity,  and  Eqs.  (20)  and  (21)  coincide  with  the 
hydrodynamic  equations  for  shallow  water  (see,  for  example,  2f).  This  means  that  the  2D-eleetron  fluid  in  a 
Ballistic  FET  should  behave  like  shallow  water.  In  this  hydrodynamic  analogy,  v  corresponds  to  the  fluid 
velocity,  and  qU/m  corresponds  to  gh,  where  h  is  the  shallow  water  level  and  g  is  the  free  fall  acceleration. 

Hence,  phenomena  similar  to  wave  and  soliton  propagation,  hydraulic  jump,  and  the  "choking"  effect 28 
should  take  place  in  a  ballistic  hydrodynamic  electron  fluid.  Effects  of  collisions,  surface  scattering,  changes  in 
the  channel  cross  section,  and  others  may  be  also  understood  using  this  analogy. 

In  the  opposite  limited  case  of  a  collision  dominated  transport,  Eq.  (20)  becomes 

qr m  dU 

v  =  ~—  -  =mE,  (20a) 

m  ox 

which  is  usually  replaced  by  a  field-dependent  velocity  function 

v  =  v(£),  (23) 

which  accounts  for  the  electron  velocity  saturation  in  high  electric  fields  and/or  for  a  region  of  negative 
differential  mobility.  Examples  of  empirical  velocity-field  dependences  are 

v  AE 

^\  +  {mE/vs)2  ^4) 

for  electrons  in  silicon,  and 


'  t  TT^ITiyj 

for  electrons  in  compound  semiconductor  materials  with  a  negative  differential  mobility  region,  such  as  GaAs. 
Here  Es  =  v5/p  is  the  velocity  saturation  field,  v5  is  the  electron  saturation  velocity,  and  A  and  t  are  fitting 
parameters.  The  analytical  model  of  a  FET,  which  accounts  for  the  saturation  of  the  electron  velocity  and  for  the 
source  and  drain  series  resistances  yields  the  following  expressions  for  the  device  drain  current 


where 


9 


(26) 


goo=Cif4nj(Ugs-UT) 

_ ScJUy  - Ut J 

( U  -U  \ 

1  +  gc>uA  + ,  1  +  2  gchoRs  +  -g  f 

I  UL  ) 


(27) 

(28) 
(29) 


Ul  -  ESL,  and  Rs  and  Rj  are  source  and  drain  series  resistances,  respectively. 

In  the  collision  dominated  regime,  operating  regime,  the  upper  frequency  of  operation  of  a  Field  Effect 
Transistor  (FET)  is  limited  by  the  electron  transit  time,  tlr.  The  transistor  cutoff  frequency,  fT  =  1/(2 ntlr),  see  for 
example,  .  However,  plasma  effects  become  important  in  modem,  short  channel  field  effect  transistors,  where 


the  sheet  carrier  density  is  very  high.  These  effects  should  allow  us  to  use  FETs  at  much  higher  frequencies  (in 
the  terahertz  range  for  deep  submicron  devices). 

A  explained  above,  plasma  waves  with  a  linear  dispersion  law,  ©  =  sk.  The  velocity  of  the  plasma  waves,  5, 
is  typically  on  the  order  of  10s  cm/s,  which  is  much  larger  than  the  drift  velocity  of  the  2D  electrons  in  the  FET 
channel.  This  is  why  the  propagation  of  plasma  waves  can  be  used  for  new  regimes  of  FET  operation,  with  a 
much  higher  frequency  than  for  conventional,  transit-time  limited  regimes.  Under  certain  conditions,  plasma 
oscillations  can  be  excited  in  a  FET  by  a  dc  current,  and  the  FET  can  be  used  as  an  oscillator  operating  in  the 
terahertz  range.  Nonlinear  properties  of  the  plasma  waves  can  be  utilized  for  terahertz  detectors,  broad  band 
detectors,  mixers,  and  frequency  multipliers. 30 

Devices  operating  at  terahertz  frequencies  should  be  considered  to  be  inseparable  from  the  circuit. 
Therefore,  issues  related  to  coupling  plasma  wave  to  electromagnetic  radiation  conversion,  to  antenna  structures 
for  electromagnetic  radiation,  and  to  device  integration  with  submillimeter  circuits  will  have  to  be  addressed  for 
practical  implementation  of  plasma  wave  electronic  devices.  Since  the  plasma  wave  velocity  is  much  smaller 
than  the  light  velocity,  and  the  device  dimensions  are  much  smaller  than  the  electromagnetic  wave  length 
corresponding  to  the  plasma  frequency,  antenna  structures  —  which  are  much  larger  than  typical  devices  --  are 
needed  for  coupling  plasma  and  electromagnetic  waves.  These  issues  and  the  antenna  and  circuit  design  will  be 
similar  to  those  currently  investigated  for  deep  submicron  Schottky  diodes  operating  in  the  terahertz  range. 3i 


4.Plasma  wave  instability 

Dyakonov  and  Shur  showed  that  the  surface  plasma  waves  propagating  in  gated  2DEG)  might  grow  in  a 
ballistic  device  due  to  reflections  from  the  device  boundaries. 16 

The  linearized  system  of  Eqs.  (20)  (in  the  limit  of  r->oo)  and  (22)  predicts  the  dispersion  law,  k=  ±afs, 
which  is  the  same  as  for  the  shallow  water  waves.  The  wave  velocity  5  =  (q LJ0/ tn) ^  If  the  electrons  move 
with  a  drift  velocity  v0,  the  waves  are  carried  along  by  the  flow,  and  the  dispersion  relation  becomes  k  =  ox(v0  ± 
s).  (The  drift  velocity,  v0,  corresponds  to  the  current  per  unit  width  j  =  qnsva  =  CU0vJq.) 

Dyakonov  and  Shur  assumed  the  following  boundary  conditions  at  the  source  and  drain  side  the  device 
channel:  constant  gate-to-source  voltage,  Ugs,  at  the  source  side  of  the  channel  (x  =  0)  and  constant  drain  current 
at  the  drain  side  of  the  channel  (x  =  L).  The  realization  of  these  boundary  conditions  is  straightforward  at  low 
frequencies,  i.e.  one  applies  voltage  between  the  gate  and  source  using  a  voltage  source,  Ugs,  and  connects  a 
current  source,  /<*,  to  the  drain.  At  high  frequencies,  the  boundary  conditions,  such  as  short  and  open  circuits, 
cannot  be  easily  implemented.  Also,  the  displacement  current  between  the  gate  and  the  channel  plays  an 
important  a  role.  However,  all  microwave,  millimeter,  and  submillimeter  field  effect  transistors  designed  in 
such  a  way  that  the  gate-to-source  parasitic  and  fringing  capacitances  (which  add  to  the  input  capacitance)  is 
greater  or  even  much  greater  than  the  gate-to-drain  parasitic  and  fringing  capacitances  (which  contribute  to  the 
Miller  capacitance).  The  intrinsic  gate-to-source  and  gate-to-drain  capacitances  are  equal  at  zero  drain  bias. 
However,  the  intrinsic  gate-to-source  capacitance  increases  with  an  increase  in  the  drain  bias  and  the  intrinsic 
gate-to-drain  capacitance  decreases  with  an  increase  in  the  drain  bias  (see  Fig.  5), 111 

Therefore,  the  boundary  conditions  at  the  source  are  closer  to  the  short  circuit  conditions  at  high  frequencies 
and  the  boundary  conditions  at  the  drain  are  closer  to  the  open  circuit.  This  asymmetry  increases  with  an 
increase  in  the  drain  bias  as  shown  in  Fig.  5.  We  let  v  =  v0  +  viexp(-iraf),  U=  Ua+  Uxexp(-icot),  linearize  Eqs. 
(20),  (22)  with  respect  to  Vi  and  Uu  and  use  the  boundary  conditions  U0)  =  0  and  Aj'(L)  =  0  (i,  e.,  U0v\(L)  + 
v0Ui(L)  =  0)  as  discussed  above.  We  then  seek  the  solution  as  the  sum  of  two  waves  propagating  from  the 
source  to  the  drain  and  from  the  drain  to  the  source,  respectfully,  with  the  wave  vectors  k\  and  ko: 


II  The  analogy  between  plasma  waves  in  a  FET  and  water  waves  goes  even  further.  While  for  kd  «  |  the  plasma  waves  with  O)  ~  k  arc  simrlar  to  the  shallow  water  waves,  in  the  opposite  limiting  case  kd 
»  1 .  the  plasma  waves  have  the  same  dispersion  law  Ai/2asthe  gravitational  waves  in  deep  water. 

1,1  We  should  notice  that  capacitances  shown  in  Fig.  5  are  calculated  assuming  the  collision  dominated  transport  in  a  FET.  However, 
the  qualitative  conclusions  should  apply  to  the  ballistic  regime  as  well. 


where 


v,  =  /4exp(i£,x)+  5exp(/£2x) 
Ux  =  C  expf/^x)  +  D  exp  (ik2  x) 
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Fig.  5.  Gate-to-source  and  gate-to-drain  capacitances  for  zero  and  -0.3  gate  biases:  Cgs  and 
CoD,ve rsus  drain  bias.  Parameters  used  in  the  calculation:  threshold  voltage,  UT=-  0.6  V;  gate 
width  and  length  W=50  mm  and  1=0.15  mm,  respectively;  AlGaAs  thickness,  d,  =  25  nm;  and 
maximum  density  of  the  two-dimensional  electron  gas,  nsmax  =2.5x1 016  m'2. 32 

This  procedure  leads  to  the  following  expressions  for  the  real  and  imaginary  parts  of  (o  =  (o'+io)" : 


s2  -v2 


where  n  is  an  odd  integer  for  jv„j  <  s  and  an  even  integer  for  |vc|  >  s.  Eq.  (35)  shows  that  for  positive  va,  the 

steady  flow  is  unstable  if  v0<s  and  stable  if  v0>s.  For  vjs  «  1 ,  co"=  vJL  which  is  the  inverse  electron 
transit  time. 

The  boundary'  conditions  play  a  very  important  role  in  determining  the  growth  of  this  plasma  instability. 
Different  boundary  conditions  might  be  interpreted  in  terms  of  effective  impedances  connected  between  the 
source  and  gate  and  the  gate  and  drain,  respectively  (see  Fig.  6): 

z  -JIM.  ™ 

-  j,(oyv  <36) 

where  j\  =  qnsXva  +  qnsov]  and  W  is  the  gate  width. 
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Fig.6.  Boundary  conditions  for  plasma  waves  in  FET  channel. 


Our  analysis  showed  that  the  instability  might  take  place  as  long  as  Jzgs| « |ZL>. j .  A  more  detailed  analysis  of 

different  boundary  conditions  for  this  system  was  given  by  Crowne. 33 

The  above  analysis  is  based  on  the  solution  of  hydrodynamic  equations,  which  are  valid  when  electron- 
electron  collisions  are  very  frequent.  Dmitriev  et  al 31  presented  the  solution  of  the  Boltzmann  equation  for  a 
relatively  low  density  gated  two-dimensional  electron  gas,  where  electron-electron  collisions  are  not  significant. 
This  solution  revealed  that  the  plasma  waves  have  the  same  dispersion  law  as  for  a  high  electron  sheet  density. 
In  both  cases,  the  plasma  waves  become  instable  in  a  short  channel  field  effect  transistor  with  asymmetric 
boundary  conditions  at  small  channel  currents.  They  also  investigated  the  role  of  the  Landau  damping.  The 
Landau  damping  is  the  damping  of  a  plasma  wave  caused  by  a  transfer  of  energy  from  the  wave  to  particles, 
whose  velocity  equals  the  phase  velocity  of  the  wave.  In  this  system,  the  electron  velocity  is  vp  (or  vth  in  a  non¬ 
degenerate  case)  and  the  plasma  wave  velocity  is  s.  Dmitriev  et  al  showed  that  the  Landau  damping  is  small 
when  s  »  vp,  which  is  quite  understandable.  The  inequality  s  »  vp  can  be  presented  as 

d»rj 2  (38) 


where  rs  = 


Ana  ft2 
mq2 


is  the  Bohr  radius.  For  GaAs,  e  =  1.1 4x1  O'10  F/m,  m  =  0.067  mQ,  where  m0  is  the  free 


electron  mass,  and  rs  =  10  nm.  Typical  values  of  d  are  on  the  order  of  30  nm,  so  that  Eq.  (38)  is  valid,  even 
though  this  estimate  shows  that  the  Landau  damping  might  be  important  in  devices  with  a  smaller  separation 
between  the  gate  and  the  channel. 


5.  Instability  conditions 

Three  decay  mechanisms  oppose  plasma  wave  growth  caused  by  the  instability  described  in  Section  4:  electron 
scattering  by  phonons  or  impurities,  the  effect  of  the  finite  electron  transit  time  in  the  channel  leading  to  the 
“ballistic  mobility”  given  by  Eq.  (2),  and  internal  friction  caused  by  the  viscosity  of  the  electron  fluid. 

The  electron  scattering  can  be  accounted  for  by  retaining  the  term  -v/xm  in  the  right-hand  side  of  Eq.  (20). 
This  leads  to  the  addition  of  the  -  1/(2 1 m)  term  to  the  wave  increment.  Hence,  the  wave  grows  only  if  the 
number  of  scattering  events  during  the  transit  time  is  small. 

The  viscosity,  v,  of  the  electron  fluid  causes  an  additional  damping  with  the  decrement  of  vk~  where  k  is  the 
wave  vector.  Hence,  the  viscosity  is  especially  effective  in  damping  higher  order  modes.  Comparing  d'  with 
vk~  for  the  first  mode,  we  find  that  the  effect  of  the  viscosity  is  small  when  the  Reynolds  number  Re  =  Lv0  N  is 
much  greater  than  unity.  In  a  highly  non-ideal  electron  gas,  where  the  Bohr  energy,  thermal  energy,  and  Fermi 
energy  are  of  the  same  order  (which  roughly  corresponds  to  the  surface  electron  concentration  of  1012  cm'2  at  77 
K  in  GaAs),  the  viscosity  of  the  electron  fluid  is  on  the  order  of  h/  m ,  which  is  approximately  1 5  cm2/s 
(comparable  to  that  of  castor  oil  or  glycerin  at  room  temperature).  The  Reynolds  number  may  be  estimated  as 
Re  =  mv0L/h~  12  for  v0  =  107  cm/s  and  L  =  0.2  pm  (see  Fig.  7). 

For  a  sample  with  L  =  0.2  pm  at  77  K,  assuming  xm  ~  1 0'1 1  s,  the  increment  vJL  exceeds  the  decrement 
l/(2xw)  caused  by  the  collisions  when  vc  >  106  cm/s.  For  the  same  sample,  the  decrement  caused  by  viscosity, 
v(2tc/Z)2/16,  is  smaller  then  the  increment  vJL  when  v  =  ;r2v/(4l)  ~  1.8xl06  cm/s.  Hence,  the  threshold 
velocity  for  the  instability  is  well  below  the  peak  velocity  in  GaAs  (~2xl07  cm/s). 


Fig.  7.  Reynolds  number  versus  gate  length  for  different  velocity  values  in  FET  channel 

The  condition  vJL<M(2xm)  corresponds  to  the  requirement  of  a  ballistic  transport  in  a  transistor  channel.  In 
this  regime, 

v0  =  MbaiUds/L  (39) 

where  Uds  is  the  intrinsic  drain-to-source  voltage  drop  across  the  channel  (excluding  the  voltage  drop  across  the 
source  and  drain  series  resistances)  and  ai  is  given  by  Eq.  (2)  for  a  non-degenerate  case.  Hence,  the  inequality 
vJL<  l/(2x  m)  can  be  re-written  as 


ud,>  Ucr  = - —  ^  (40) 

M 

Fig.  8  shows  the  values  of  Ucr  and  vocr  =  L/(2rm)  for  GaAs-based  HEMTs  versus  the  gate  length  at  300  K  and  77 
K,  As  seen  from  these  figures,  the  instability  condition  can  be  achieved  in  deep  submicron  (less  than  0.1  micron 
or  so)  devices  at  room  temperature  and  in  one  to  two  micron  size  (or  shorter)  devices  at  77  K. 

The  first  measurement  results  by  Cheremisin  point  out  that  such  instability  might  indeed  occur  in  0.1  micron 

CaAs-based  HEMTs.  3s  However,  more  experimental  studies  are  needed  to  prove  the  existence  of  this 
instability. 

Recently  Ryzhii  and  Shur  proposed  to  combine  this  approach  with  using  resonant  tunneling  structure 
inserted  between  the  gate  and  a  channel  of  HEMT  or  a  related  NERFET-type  device  ^  (see  Fig.  9),  Their  theory 
predicts  a  very  large  increase  of  the  plasma  wave  instability  increment.  Ryzhii  et  aL*  also  proposed  to  use  a 
similar  structure  to  enhance  the  plasma  wave  detector  sensitivity  (see  Section  6). 
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Fig.  8.  Ucr  and  vocr  =  L/(2tp)  for  GaAs-based  HEMTs  versus  the  gate  length  at  300  K  and  77  K 


Fig.  9.  Resonant  tunneling  structures  for  enhancement  of  plasma  wave  instability. 


For  a  FET  with  a  resonant  tunneling  structure  between  the  gate  and  the  channel,  Eq.  (21)  should  be  re-written  as 

<?(n,v)  «  .  ( eRT  1 

— L+--  = — an 


where  re  is  the  product  of  the  try-to-eseape  frequency  and  the  maximum  transmission,  A(z)  =  — iy  is  the 

1  +  z~ 

resonance  tunneling  form-factor  ,  ert  =  Ert  -qaU/2  and  F  are  the  energy  and  position  of  the  resonant  tunneling 
level,  and  a  ~  1  is  a  geometrical  factor. 

The  boundary  conditions  for  this  structure  were  chosen  as 

U(0)  =  U(l)-0  (42) 

As  was  shown  in  Reference  j3,  the  plasma  wave  frequencies  in  such  system  are  given  by 

(43) 

where  k„  =  nn/L,  and  n  =  1 , 2,  3  ... 


Fig.  10  shows  the  increment  of  the  instability  calculated  in  Reference  for  0.5  micron  gate  devices.  As 
seen,  this  instability  should  take  place  in  submicron  GaAs-based  FIEMTs  with  a  resonant  tunneling  structure  at 
77  K,  when  the  electron  mobility  is  large  enough. 

Ryzhii  mid  Shur  also  suggested  that  a  forward  gate  current  in  FETs  can  enhance  the  plasma  instability 
increment.  Electrons  constituting  this  current  experience  the  time  delay  that  typically  might  be  on  the  order  of 
d/veff,  where  vejf  is  the  velocity  of  the  electrons  crossing  from  the  channel  into  the  gate.  Such  delay  might  lead  to 
a  dynamic  negative  conductance  at  plasma  frequencies,  which  should  result  in  the  excitation  of  plasma  waves  in 
the  transistor  channel. 


Once  the  electron  velocity  exceeds  the  threshold,  the  plasma  waves  grow.  Since  no  other  steady  states  exist 
for  v0  <  s,  this  growth  should  lead  to  oscillations  for  which  the  plasma  wave  amplitude  is  limited  by  non¬ 
linearity.  Nonlinear  plasma  oscillations  in  FETs  have  been  studied  in  References4*3,4*’  . 

Let  us  now  discuss  possible  applications  of  this  instability.  The  plasma  oscillations  result  in  a  periodic 
variation  of  the  channel  charge  and  the  mirror  image  charge  in  the  gate  contact,  i.  e.  to  the  periodic  variation  of 
the  dipole  moment.  This  variation  should  lead  to  electromagnetic  radiation.  The  device  length  is  much  smaller 
than  the  wavelength  of  the  electromagnetic  radiation,  XR,  at  the  plasma  wave  frequency.  (The  transverse 
dimension,  W,  may  be  made  comparable  to  XR).  Hence,  the  Ballistic  FET  operates  as  a  point  or  linear  source  of 
electromagnetic  radiation.  Many  such  devices  can  be  placed  into  a  quasi-optical  array  for  power  combining. 
This  device  de-couples  the  operating  frequency  (which  can  be  tuned  in  a  wide  frequency  range  by  varying  IJ0 ) 
from  the  electron  transit  time  limitation.  The  maximum  modulation  frequency  is  still  limited  by  the  transit  time. 


Fig.  10.  Increment  of  fundamental  mode  versus  channel  mobility  for  resonant  tunneling  HETM 
structures. 33 


6.  Detectors  and  mixers  of  terahertz  radiation 

The  excitation  of  surface  glasma  waves  can  be  also  used  for  detection,  multiplication,  and  mixing  of  terahertz 
radiation.  Weikle  et  al.  and  Lu  et  al.  43  reported  on  the  detection  of  the  microwave  and  terahertz  radiations, 
respectively,  using  AlGaAs/GaAs  HEMTs  Fig.  1 1  shows  a  schematic  diagram  and  an  equivalent  circuit  of  a 
HEMT  subjected  to  terahertz  radiation. 


IV  We  should  notice  that  for  M  >  S,  in  addition  to  the  uniform  flow,  there  is  another  steady  state  solution  corresponding  to  a  hydraulic  jump  in  shallow  water  16 
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Fig.  1 1 .  Schematic  geometry  of  FET  operating  in  detector  mode  for  induced  ac  voltage  (a)  and 
induced  ac  current  (b). 15  ©  IEEE,  1998 


The  incoming  electromagnetic  radiation  induces  an  ac  voltage  at  the  source  side  of  the  channel.  The  drain  side 
of  the  channel  is  an  open  circuit.  (A  coupling  with  a  terahertz  radiation  might  be  enhanced  using  an  appropriate 
antenna  structure.)  This  ac  voltage  excites  the  plasma  waves.  As  mentioned  above,  the  dc  gate-to-source 
voltage  determines  the  velocity,  s,  of  these  waves.  Because  of  the  nonlinear  properties  of  the  electron  fluid  and 
the  asymmetry  in  the  boundary  conditions,  a  FET  biased  only  by  the  gate-to-source  voltage  and  subjected  to 
electromagnetic  radiation  develops  a  constant  drain-to-source  voltage,  which  has  a  resonant  dependence  on  the 
radiation  frequency  with  maxima  at  the  plasma  oscillation  frequency  and  its  odd  harmonics.  The  induced  drain 
voltage,  A  U,  is  given  by  15 


where 
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and  ka  and  k0  are  the  real  and  imaginary  parts  of  the  wave  vector,  k0: 
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Here  where  n  l,  3,  5,  7...,  U0  is  the  gate-to-source  gate  voltage  swing,  and  Ua  is  the  amplitude  of  the  ac 
source-to-gate  voltage  induced  by  the  incoming  radiation.  The  half  width  of  the  resonance  curve  is  determined 


by  the  damping  of  the  plasma  oscillations  caused  by  the  electron  momentum  relaxation  and/or  the  electron  fluid 
viscosity.  Thus,  the  FET  acts  as  a  tunable  resonance  quadratic  detector  of  electromagnetic  radiation. 

Lu  et  al  reported  on  the  implementation  of  a  terahertz  detector  utilizing  2D  electronic  fluid  in  a  High 
Electron  Mobility  Transistor  (HEMT). 24  To  our  knowledge,  this  is  the  first  plasma  wave  terahertz  detector  ever 
demonstrated.  The  detector  was  fabricated  using  a  Fujitsu  FHR20X  HEMT  mounted  on  a  quartz  substrate.  The 
device  operated  at  a  frequency  of  2.5  THz,  which  is  about  30  times  higher  than  the  transistor  cutoff  frequency, 
A  CC>2-pumped  far-infrared  gas  laser  served  as  a  source  of  terahertz  radiation.  The  laser  beam  was  chopped  and 
focused  on  the  sample  with  the  electric  field  polarization  oriented  in  the  drain-to-source  direction.  The  drain 
was  open,  thus  the  drain  current  is  zero.  The  gate  current  is  in  the  order  of  nano-ampere  for  the  measurement 
range  of  gate  bias.  In  agreement  with  the  predictions  of  our  terahertz  detector  theory  ,  the  radiation  induced  a 
negative  DC  drain-to-source  bias  proportional  to  the  radiation  intensity.  This  voltage  was  measured  using  a 
lock-in  amplifier.  The  measured  dependencies  of  the  detector  responsivity  on  the  gate  bias  are  in  good 
agreement  with  the  gate  bias  dependence  of  the  normalized  responsivity  predicted  by  the  detector  theory  15  (see 
Fig.  12).  The  responsivity  increases  at  smaller  gate  voltage  swings  as  predicted  by  the  theory.  The 
dimensionless  responsivity  function  displays  a  rather  broad  resonant  peak  corresponding  to  the  third  harmonic 
of  the  surface  plasma  frequency. 

More  recently,  Knap  et  al  44  demonstrated  the  resonant  detection  of  sub-terahertz  radiation  by  the  two- 
dimensional  electron  plasma  confined  in  a  submicron  gate  GaAs/AlGaAs  field-effect  transistors.  Their  results 
confirmed  that  the  critical  parameter  governing  the  sensitivity  of  the  resonant  detection  is  cojm,  where  co  is  the 
radiation  frequency  and  vm,  is  the  momentum  scattering  time.  By  lowering  the  temperature  and  hence 
increasing  zm  and  increasing  the  detection  frequency  co,  they  reached  the  mm  ~  1  condition  and  observed  a 
resonant  detection  of  600  GHz  radiation  in  the  0.15  pm  gate  length  GaAs  field-effect  transistor  (see  Fig.  13). 
The  evolution  of  the  observed  photoresponse  signal  with  temperature  and  frequency  is  well  reproduced  within 
the  framework  of  a  theoretical  model  developed  earlier. 

Let  us  now  consider  an  HFET  detector  with  a  longer  channel,  such  that  sxJL  «  1.  In  this  case,  the  plasma 
waves  are  excited  by  the  incoming  radiation  near  the  source,  provided  that  cox  »  1  but  they  decay  before  ever 
reaching  the  drain  side  of  the  channel. 
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Fig.  12.  Measured  (symbols)  response  versus  incident  power  at  zero  gate  bias  at  2.5  THz  and  the 
laser  power  of  5.7  mW.  The  solid  line  is  normalized  for  comparison.  The  peak  corresponds  to  the 
third  harmonic  of  the  surface  plasma  frequency,  e>3.’32  ©  IEEE,  1 998 
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Fig.  13  a)  Experimental  drain  response,  Uds,  for  different  temperatures.  Inset  shows  the  resonant 
signal  AUds  at  8  K  obtained  after  subtraction  of  the  1/Uo  like  background,  b)  Resonant  plasma 
frequency  as  a  function  of  the  gate  voltage.  The  vertical  arrow  shows  the  cross  point  with 
600GHz  horizontal  line  —  the  voltage  of  the  expected  resonant  detection. 44 
The  analysis  of  this  case,  based  on  the  solution  of  Eqs.  (44)  -  (48),  shows  that  the  dc  voltage  will  still  be 
developed  between  the  drain  and  source,  and  the  device  will  operate  as  a  broad  band  detector  of 
electromagnetic  radiation  with  the  output  signal 


(49) 


As  can  be  seen  from  Eq.  (49),  a  long  HFET  should  act  as  a  broad  band  detector  of  electromagnetic  radiation. 
The  highest  frequency  of  detection  is  on  the  order  of  s/d  where  d  is  the  gate-to-source  spacing,  since  the  gradual 
channel  approximation  is  not  valid  for  the  plasma  wavelengths  shorter  than  d.  For  s  ~  108  cm/s  and  d~  100  A, 
this  corresponds  to  1 00  THz. 

Weikle  et  al  fabricated  a  prototype  non-resonant  detector  (operating  in  the  microwave  range)  using  an 
AlGaAs/GaAs  0.15  micron  gate  HEMT.  4~  The  measured  dependencies  of  the  detector  responsivity  on  the  gate 
bias  and  frequency  were  in  good  agreement  with  our  theory.  AlGaN/GaN  HFETs  exhibited  a  non-resonant 
detection  at  frequencies  much  higher  than  the  HEMT  cutoff  frequency.  The  measured  detector  responsivity  was 
in  good  agreement  with  the  theory  and  reaches  300  V/W,  which  is  comparable  to  the  responsivity  of  Schottky 
diode  detectors  at  these  frequencies.  This  value  can  be  greatly  improved  by  using  a  proper  antenna  structure. 

More  recently.  Knap  et  al  4i  reported  on  an  experimental  and  theoretical  study  of  nonresonant  detection  of 
sub-terahertz  radiation  in  GaAs/AlGaAs  and  GaN/AlGaN  heterostructure  field  effect  transistors  (see  Fig.  14), 
The  experiments  were  performed  in  a  wide  range  of  temperatures  (8-300K)  and  for  frequencies  ranging  from 
100GHz  to  600GHz.  The  photoresponse  measured  as  a  function  of  the  gate  voltage  exhibited  a  maximum  near 
the  threshold  voltage.  The  results  were  interpreted  using  a  new  theoretical  model  that  shows  that  the  maximum 
in  photoresponse  and  showed  explained  by  the  combined  effect  of  exponential  decrease  of  the  electron  density 
and  the  gate  leakage  current. 


Gate  Voltage  Ug  (V) 


Fig.  14,  Response  at  600  GHz,  measured  dots,  calculated  solid  lines  (a)  and  drain  current  versus 
gate  voltage  Ug  (b)  measured  in  three  experiments  Tl,  12,  T3.  Results  marked  as  T1  corresponds 
to  the  transistor  with  the  threshold  voltage  Uth=~Q.55V  measured  at  300K.  Results  marked  as  T2 
corresponds  to  the  transistor  with  Uth  =-0.42V.  Results  marked  as  T3  correspond  to  the  same 
device  but  measured  at  temperature  of  10  K,  for  which  the  threshold  voltage  was  lower  (U,h  =- 
0.22V).45 

These  detectors  operated  at  frequencies  well  above  the  cutoff  frequency.  Knap  et  al  obtained  similar 
experimental  results  for  the  broadband  sub-terahertz  detectors  fabricated  using  AlGaN/GaN  HFETs.  Their 
results  showed  that  the  2D  electron  plasma  effects  are  indeed  universal,  not  dependent  upon  a  particular 
material  system  (see  also  46) 

7.  Terahertz  photomixing  using  resonant  excitation  of  plasma  oscillations 

Alternative  approaches  to  generate  THz  radiation  are  associated  with  optical  techniques,  that  use  a  coherent 
output  at  the  difference  frequency  (equal  to  the  difference  between  the  frequencies  of  radiation  emitted  by  two 
lasers)  or  a  response  of  photoconductive  structures  to  femtosecond  optical  pulses.42"5*  Fast  quantum  well 
infrared  photodetectors  (QWIPs)  utilizing  intersubband  transitions 52  can  also  be  used  for  the  generation  of  THz 
radiation  by  mixing  infrared  laser  beams.  Indeed,  as  shown  theoretically,  QWIPs  can  exhibit  a  marked  response 
to  infrared  signals  in_the  THz  range  of  modulation  frequencies  if  the  electron  transit  time  is  short  enough  (as  can 
be  in  single  QWIPs  5 ,  )  or  if  electrons  reveal  a  pronounced  velocity  overshoot  after  their  photoexcitation  from 
QWs.  The  THz  signals  produced  by  coherent  plasma  oscillations  of  the  photogenerated  carriers  have  been 
observed  in  p-i-n  structures  by  Sha  et  al.S7  and  Kersting  et  al.  58  However,  due  to  a  strong  damping  of  the 
plasma  oscillations  excited  by  short  optical  pulses,  only  few-cycle  THz  signals  have  been  observed. 

Recently,  Ryzhii  et  al  predicted  that  modulated  infrared  radiation  can  cause  the  resonant  excitation  of 
plasma  oscillations  in  quantum  well  diode  and  transistor  structures  (see  Fig.  15  and  16).  This  effect  provides  a 
new  mechanism  for  the  generation  of  tunable  terahertz  radiation  using  photomixing  of  infrared  signals.  They 
developed  a  device  model  for  a  Quantum  Well  Photo  Mixer  (QWPM)  and  calculated  its  high-frequency 
performance.  It  was  shown  that  the  proposed  device  can  significantly  surpass  photomixers  utilizing  standard 
quantum  well  infrared  photodetectors. 


Fig.  15.  Schematic  view  of  diode-type  (a)  and  transistor-type  (b)  QWPM  structures  without 
doping  of  the  gate  barrier  and  a  band  diagram  of  transistor-type  QWPM  (c).  Arrows  indicate 
electron  trajectories  (bound-to-continuum  photoexcitation  followed  by  transport  over  collector 
barrier). 6 

Ryzhii  et  al  also  proposed  a  QWPM  based  on  QW  diode  and  transistor  structures  utilizing  bound-to- 
continuum  electron  transitions,  which  can  be  used  for  photomixing,  and  evaluated  the  device  operation  using  an 
analytical  model.  It  was  shown  that  the  QWPM  with  high  electron  mobility  in  the  QW  channel  can  exhibit  a 
resonant  response  to  modulated  infrared  radiation  when  the  modulation  frequency  is  close  to  one  of  the  plasma 


frequencies  which  can  be  in  the  THz  range.  High  values  of  the  responsivity  in  the  THz  range  of  signal 
frequencies  indicated  that  the  proposed  QWPM  and  arrays  of  such  devices  can  be  used  for  efficient  generation 
of  THz  radiation  using  photomixing  of  infrared  signals. 
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FIG.  16.  Characteristic  photoelectric  gain  vs  signal  frequency  for  QWPMs  with  different  electron 
mobilities  (different  Q  factors)  in  quantum  well  channel. 6 

8.  Conclusions. 

Plasma  wave  excitation  in  submicron  field  effect  transistors  and  related  device  structures  should  allow  us  to 
develop  a  new  generation  of  solid-state  terahertz  tunable  terahertz  devices  that  will  find  numerous  applications 
in  industry,  defense,  and  biotechnology.  Recent  experimental  results  demonstrated  both  resonant  and  non¬ 
resonant  detection  of  terahertz  radiation  by  plasma  waves.  These  results  are  quite  encouraging  and  indicate  that 
proposed  improved  structures  will  be  viable.  New  improved  structures  will  use  an  enhancement  of  the  plasma 
wave  growth  due  to  static  or  dynamic  negative  differential  conductivity  induced  by  the  negative  differential 
conductivity  caused  by  the  gate  current. 
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